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Abstract. We study the fluid/solid phase transition via a mean field model. Our first 
result is the derivation of a variational characterization of the entropy density, compatible 
with the infinite volume limit. We then determine the optimizing particle distributions for 
small energy density and a range of particle density, though we can only prove they are 
local, not global, minima of the entropy density. With this assumption we then prove that 
the resulting entropy density must lose its analyticity in various regimes. In particular this 
implies the existence of a phase transition between distinct heterogeneous structures at low 
energy density, and a phase transition between these structured phases and the disordered 
phase at high energy density. 



1. Introduction 

We will adapt an old random graph model of Strauss [St] to provide a mean field model 
of the fluid/solid phase transition of equilibrium matter, in particular to model how high 
density can produce the internal (crystalline) structure of solids out of a homogeneous fluid. 
The energy ground state in the model exhibits a range of different structures as the particle 
density varies, and we use this feature to show (to second order of perturbation theory) 
that there is a (solid/solid) phase transition between the structures, near the energy ground 
state, as well as a (fluid/solid) transition between the structured phases and the high energy 
disordered phase. We are forced to use a mean field model as there is still no convincing 
model of the basic phenomenon of a distinct solid phase; see |A1 |Brl IUh[ EH lARlj . We 
emphasize that ours is a toy model, in which the energy is not based on the physics of any 
real material but only has vaguely appropriate features of repulsion and attraction. 

The following is the usual route for understanding the thermodynamic phases of a noble gas 
such as argon, using classical statistical mechanics and a Lennard- Jones, rotation symmetric, 
2-body interaction potential V, depending on separation r by V{r) = cir~^^ — C2r~^, with 
ci,C2 > 0; see Figure 1. This phenomenological interaction contains the desired features 
that the (neutral) atoms repel strongly at small separation due to their electron clouds, and 
deform at intermediate separation, providing a weak attraction responsible for the 'molecular 
bonds' in the solid phase. The two basic phase transitions, the gas/liquid and fluid/solid, 
are drawn schematically in Figure 2; neither can be proven for this model, and indeed it has 
yet to be proven even that the energy ground state is crystalline. (The appropriate crystal 
has been proven to be the unique ground state in 1 dimension jGR] and to be a ground state 
in 2 dimensions jT].) 
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Figure 1. The Lennard- Jones interaction potential 

Other, less-reahstic models supply a convincing theoretical understanding for a gas/liquid 
transition for such materials. This was first achieved by the well known mean field model 
of van der Waals, in which the attraction between particles is simplified by ignoring the 
separation. See |LMP] for recent history of the model and improvements beyond mean field. 
These models reproduce the basic features of the transition as first order, with coexisting 
disordered phases of different energy and particle density, ending in a critical point; see 
Figure [21 

We note that not just for noble gases, but also for most materials, the gas/liquid transition 
is associated with an attractive force between the molecules. For argon the Lennard- Jones 
potential is reasonable, but for molecules which utilize other types of bonds in the solid phase 
a more complicated modeling would be desired. In contrast, the fluid/solid transition, at 
least at high pressure, is normally associated with the repulsion of the electron clouds of the 
molecules, independently of the detailed interaction associated with the outer electrons. The 
hard sphere model [Low], in which the only interaction is a simple hard core — two particles 
separated by less than the hard core distance produce infinite energy — is then appropriate 
quite generally, at least for roundish molecules. Simulation of the hard sphere model seems 
to show only a densely packed face centered cubic crystal ( [BFMHl IW]): the less-dense 
crystals which appear at lower pressure in some materials presumably are in part due to 
other features of the interaction, as in the simple cubic crystals of ionic solids such as table 
salt. Of course the hard sphere model does not show a gas/liquid transition since it does 
not include an attractive force. But the model does exhibit the basic phenomenon whereby 
(crystalline) structure is produced at high density. 

Unfortunately our knowledge of the hard sphere model [Lowj . and more generally the 
theoretical understanding of the creation of crystalline structure at high density, is based 
almost completely on computer simulation |Br[ lUhj . We will try to remedy this via an 
analogue of the van der Waals result, a mean field analysis but now for the fluid/solid 
transition. Following |PN[ ICD[ lARlj IAR2j we adapt a model of Strauss [St], with particles 
represented by edges in a large abstract graph and the (total) energy of a graph being the 
number of triangular subgraphs, providing a repulsion. Obviously this energy is not based 
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Figure 2. A schematic pressure/temperature phase diagram 



on the physics of real interacting molecules. The features we use are the 'point' hard core, 
associated with the restriction that between given vertices no multiple edges are allowed, 
together with the soft repulsion associated with each triangle, a 3-body interaction. (As in 
the hard sphere model there is no attraction in our model; see however |PN[ ICD[ IRY] and 
references therein for a similar approach using attraction to model a gas/liquid transition.) 
Our argument depends on several recent results on the asymptotics of large graphs. We 
begin with notation; see [LSH |LS2l IBCLSVl lLS3l ICDl ICV] for background. 

Consider the set of simple graphs G with set V{G) of (labelled) vertices, edge set 
E{G) and triangle set T{G), where the cardinality |V^(G')| = n. ('Simple' means the edges 
are undirected and there are no multiple edges or loops.) Let Z^f be the microcanonical 
partition function, the number of such graphs such that: 



e(G)^ 



e {e-6,e + 6) and t{G) = 



\T{G)\ 
(3) 



e {t-6,t + 6). 



Graphs in U„>i(ji'" are known to have edge and triangle densities, (e, t), dense in the region R 
in the e,t-plane bounded by three curves, Ci : (e,e^/^), < e < 1, the line li : (e,0), < 
e < 1/2 and a certain scalloped curve (e, /(e)), 1/2 < e < 1, lying above the curve 
(e, e(2e — 1), 1/2 < e < 1, and meeting it when e = = k/{k + 1), k > 1; see [PR] and 
references therein, and Figure 3. 

We are interested in the relative density of graphs in R, more precisely in the entropy, the 
exponential rate of growth of the number of particle configurations (graphs) as n grows, as 
follows. First consider 



(2) 



n,5 



and s 



e,t 



lim lim s 

S—>-0+ n— >oo 



n,5 
e't ■ 



(The existence of the double limit will be proven.) We will measure the growth rate by 
the entropy density t, and the main question of interest for us is the existence of phase 
transitions {i.e. lack of analyticity of Se,t) near the lower boundary of R in Figure 3. The 
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scallop loop of graph of /(e) 
t = e(2e - 1) 



(1/2,0) 
edge density e 



Figure 3. The microcanonical phase space R, outUned in sohd Unes 

lower boundary consists of the scalloped curve together with the 'first scallop', the line from 
(0,0) to (1/2,0). 

The analysis of phase transitions in traditional models with short range interactions re- 
quires the mathematical tools of the infinite volume limit. In this mean field graph setting 
appropriate tools have been developed recently under the title of 'graph limits', utilizing 
'graphons', as we sketch next. 



2. Graphons 

Consider the set W of all symmetric, measurable functions 
(3) g ■.{x,y)e [0,1]'' ^g{x,y)e [0,1]. 

Think of each axis as a continuous set of vertices of a graph. For a graph G e one 
associates 



(4) 



1 if ( \nx] , \ny] ) is an edge of G 
otherwise. 



where \y~\ denotes the smallest integer greater than or equal to y. For g E W and simple 
graph H we define 



(5) 



t{H,g)= / TT g{xi,Xj)dxi - ■ -dxi. 



where i = \V{H)\, and note that for a graph G, t{H,g'^) is the density of graph homomor- 
phisms H ^ G: 



(6) 



|hom(ij", G)\ 

\v{G)\\nH)\- 



We define an equivalence relation on W as follows: / ~ (? if and only if t{H, f) = t{H, g) for 
every simple graph H. Elements of W are called "graphons", elements of the quotient space 
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W are called "reduced graphons" , and the class containing G W is denoted g. The space 
W is compact in the metric topology with metric: 

(7) = 

where {Hj} is a countable set of simple graphs, one from each graph-equivalence class. 
Equivalent functions in W differ by a change of variables in the following sense. Let S 
be the space of measure preserving bijections a of [0,1], and for / in W and cr G S, let 
fa{x, y) = f{a{x), cr{y)). Then f g ii and only if g = for some a G S. Note that if each 
vertex of a finite graph is split into the same number of 'twins', each connected to the same 
vertices, the result stays in the same equivalence class, so for a convergent sequence g'^^ one 
may assume — )■ oo. 

The value of this formalism here is that one can use large deviations on graphs with 
independent edges [CVj to give an optimization formula for Se,f, which allows us to analyze 
Se,t near the energy ground states, the lower boundary of R in Figure 3. We next use the 
large deviations rate function for graphs with independent edges |CV] to give a variational 
characterization for the entropy density. (There is a variational characterization in [CDj of 
the analogous quantity for exponential random graph models.) 



3. A VARIATIONAL CHARACTERIZATION OF THE ENTROPY DENSITY 

Theorem 3.1. For any possible pair {e,t), Se,t = ~Hiin/((7), where the minimum is over 
all graphons g with e{g) = e and t{g) = t, where 

(^ia) = / / g{x,y)dxdy, t{g) = / / / g{x,y)g{y, z)g{z,x) dx dydz 
Jo Jo Jo Jo Jo 

and the rate function is 

(8) I{g) = I I lQ{g{x,y))dxdy, where Io{u) = l-[uln{u) + {l-u)\n{l-u)]. 

Jo Jo ^ 

Proof. We first prove that Se,t is well-defined. A priori all we know is that liminf ln(Z"f 
and limsupln(Z"f )/?T,^ exist as — t- oo. However, we will show that they both approach 
— min I{g) as 5 — ?■ 0'*'. 

We need to define a few sets. Let Us be the set of graphons g with e{g) and t{g) strictly 
within 5 of e and t, i.e. the preimage of an open square of side 26 in (e, t)-space, and let Fs 
he the preimage of the closed square. Let Us and Fs he the corresponding sets in W. Let 
\U^\ and \F^\ denote the number of graphs with n vertices whose checkerboard graphons (jl]) 
lie in Us or Fs. The large deviation principle. Theorem 2.3 of |CV] . implies that: 

In \F"'\ 

(9) limsup ' / ' < - inf I{g), 



which also equals — inf lig), and that 



In I r/" I 

(10) liminf > - inf I{g), 
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which also equals — inf I{g). This yields a chain of inequalities 
(11) 

In I I In I f/" I In I F''^ I 
— inf /(fif) < liminf < limsup < limsup < — inf I(q) < — inf I(q) 

Us n2 Fg ' ' c/^+^a 

As 5 — )■ 0"^, the limits of — inf I{g) and — inf I{g) are the same, and everything in between 

Us U^_^_^2 

is trapped. 

So far we have proven that 

(12) Se,t = - lim inf/((7). 

(5-J-0+ Us 

Next we must show that the right hand side is equal to —mml{g). By definition, we can 
find a sequence of reduced graphons gs G Us such that liml{gs) = liminf /(^f). Since W is 

<5->0 Us 

compact, these reduced graphons converge to a reduced graphon go, represented by a graphon 
go G Fo. Since / is lower- semicontinuous jC V] , I{go) < lini/(5'5), so mining) < liminf I{g). 

Fo Us 

(We write min rather than inf since Fq is compact.) However, min I{g) is at least as big as 
mil{g), since Fq C Us. Thus min/(5f) = liminf /((/). □ 

Us Fo (5-+0 C/^ 

4. Minimizing the rate function on the boundary 

From now on we will work exclusively with graphons rather than with graphs. From 
Theorem 13. all questions boil down to "minimize the rate function over such-and-such 
region" . The first region we study is the lower boundary of (e, t)-space, beginning with the 
first (flat) scallop: 

Theorem 4.1. If e < 1/2 and t = 0, then mml{g) = /o(2e)/2, and this minimum is 

Fo 

achieved at the graphon 

,^ , , \ 2e it X < - < y or y < - < x: 

(13) go{x,y) = l ^ 2 ^ ^ 2 ' 

I otherwise. 

Furthermore, any other minimizer is equivalent to go, corresponding to the same reduced 
graphon. 

Proof. Since t{g) is identically zero, the measure of the set {{x,y) E [0, l]'^\g{x,y) = 0} is at 

least 1/2. Otherwise, the graphon g{x,y) = I ' 9{ jU) > .^q^j^ have no triangles and 

1 otherwise, 

would have edge density greater than 1/2, which is impossible. So we restrict attention to 
graphons that are zero on a set of measure at least 1/2. From the convexity of /q, we know 
that the minimum of I among such graphons must be zero on a set of measure 1/2 and must 
be constant on the rest . Thus go is a minimizer, and min I{g) = I{go) = -^o(2e)/2. 

Fo 

Now suppose that g is another minimizer. Since g is zero on a set of measure 1/2 and 
is 2e on a set of measure 1/2, (7 is 1 on a set of measure 1/2, and so describes a graphon 
with edge density 1/2 and no triangles. This means that g describes a complete bipartite 
graph with the two parts having the same measure. That is, g is equivalent to the graphon 
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that equals lifx<-<yoTy<-<x and is zero everywhere else. But then g = 2eg is 
equivalent to go. □ 



, then any graph 



The situation on the curved scallops is slightly more complicated. Pick an integer i > 1. 

(The case i = 1 just gives us our first scallop.) If e G 1 — -, 1 — ^ _^ ^ 

G with edge density e and the minimum number of triangles has to take the following form 
(see |PR] for the history). Let 

There is a partition of {1, . . . ,n} into i pieces, the first £ — 1 of size [cnj and the last of 
size between [cnJ and 2 [cn\ , such that G is the complete ^-partite graph on these pieces, 
plus a number of additional edges within the last piece. ( [y\ denotes the largest integer 
greater than or equal to y.) These additional edges can take any form, as long as there are 
no triangles within the last piece. 

This means that, after possibly renumbering the vertices, the graphon for such a graph can 
be written as an uneven ixi checkerboard obtained from cutting the unit interval into pieces 
Vk = [{k — l)c, kc] ioT k < i and Ve = [(£ — l)c, 1], with the checkerboard being 1 outside the 
main diagonal, on the main diagonal except the upper right corner, and corresponding to 
a zero-triangle graph in the upper right corner. 

Limits of such graphons in the metric must take the form 

{1 X < kc < y or y < kc < X for an integer k < i; 

{k — l)c < x,y < kc for some integer k < i; 

unspecified x,y > {i — l)c, 

with 

(16) jjj g{x,y)g{y,z)g{z,x) dxdydz = 0, 

[(^-1)C,1P 

aud wth jla(.,y)d.dy = a. Minimizing /(j) on such graphons is easy, smce all but the 



[0,1]2 

upper right corner of the graphon is fixed. Applying Theorem 14.11 to that corner, we get 

Theorem 4.2. If e > 1/2 andt is the smallest value possible, then the minimum of I{g) on 

Fq is achieved by the graphon 

(17) 

{1 X < kc < y or y < kc < X for an integer k < i; 
p {i-l)c<x <[l + {i- l)c]/2 <yor{i-l)c<y<[l + {e- l)c]/2 < x; 
otherwise, 

where 

, , 4c(l - tc) 

(18) p- ^ ^ 



{\-{i-\)cf 
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is a number chosen to make 



/ / gi^iV) dx dy = e. Furthermore, any other minimizer is 
J 7[o,i]2 



equivalent to qq. The minimum value of I{g) is 

(19) I{9o) = (1 - - ^)^)' /o(p). 

5. Minimizing near the first scallop 

Now that we know the minimizer at the boundary, we perturb it to get a minimizer near 
the boundary. 

Theorem 5.1. Pick e < 1/2 and e sufficiently small. Then the graphon 

( 1 1 

/^^\ / \ I2e — e x<-<yory<-<x 

(20) g{x,y)=^ 2 ^ ^ 2 

I e otherwise, 

minimizes the rate function to second order in perturbation theory among graphons with 
e{g) = e and t{g) — — (e — e)^. For pointwise small variations Sg of g, the second 

variation in I{g) is bounded from below by - J J {5g{x,y)Y dx dy. 



[0,1 



2 



Proof. Wc first consider the first variation in I{g) for general graphons and derive the Euler- 
Lagrange equations. It is easy to check that our specific g satisfies these equations. We then 
consider the second variation in I{g). Note that the function Iq satisfies 



(21) r,{u)^\Mu)-\n{l-u)l I'M^l 



1 1 

u \ — u 



> 2. 



To find the Euler-Lagrange equations with the constraints that {e{g),t{g)) are equal to 
fixed values (cq, ^o); we use Lagrange multipliers and vary the function I{g) + Xi{e{g) — Cq) + 
^2{t{g) — to). To first order, the variation with respect to g is 



(22) SI{g) = I'Q{g{x,y))Sg{x,y)dxdy + Xi / / dg{x,y)dxdy 

Jo Jo Jo Jo 



1 /.I 



(23) +3A2 / / h{x,y)6g{x,y)dxdy, 



OJO 



where we have introduced the auxiliary function 

(24) h{x,y)= / g{x,z)g{y,z)dz. 

Jo 

Setting SI{g) equal to zero, we get 

(25) I'Mx,y)) = -Xi~3X2h{x,y). 
Our particular g{x,y) satisfies this equation with 

(26) 3A - ^0(2^-^) -TO 

^^^^ - — — ■ 
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Next we expand 5t and 51 to second order in 5g, ignoring 0{{5gf) terms. 



51^ 



J J I'o{9{x,y))59{x,y)dxdy 

+ Io{9{x,y)){6g{x,y)fdxdy 

J J (-Ai -3X2h{x,y))Sg{x,y)dxdy 

+IJJ Io{9{x,y)){Sg{x,y)fdxdy 

—Xi5e — X2St + 3A2 J J J g{x, y)5g{x, z)dg{y, z) dx dydz 

+IJJ Io{9{x,y)){Sg{x,y)fdxdy 



3A, 




9{x, y)Sg{x, z)5g{y, z)dx dydz 



(27) 



since 



+ 



\ J J Ioi9{x, y))S9{x, yfdx dy + ^JJ Io{g{x, y))Sg{x, yfdx dy, 



(28) dt = 3jj h{x,y)Sg{x,y)dxdy + 3 JJJ g{x,y)Sg{x, z)Sg{y, z)dxdydz + 0{{dgf), 



and since we are holding e{g) and t{g) fixed. We have spht the j j I"6g^ term into two 

pieces, as we will be applying different estimates to each piece. 

Since h{x, y) and I"{g) are piecewise constant, all of our integrals break down into integrals 
over different quadrants. Let Ri and R2 be the following subsets of [0, 1]^: 

(29) i?i = {x,y< l/2}U{x,y> 1/2}, R2 ^ {x < 1/2 < y}U {y < 1/2 < x}. 

/■1/2 /■! 
For each z, we define the functions fi{z) — I Sg{x, z)dx and f2{z) — I 5g{x, z)dx. The 

Jo Jl/2 



second variation in / is then 

[0,1]2 



-^jj I'M5g{x.yfdxdy+^-^ 



I"(2e — e) 

Sg{x,yfdxdy + — — // 5g{x,yfdxdy 



Ri 



R2 



+ 3A2 / dz 



e jj Sg{x,z)Sg{y,z)dxdy+{2e-e) JJ Sg{x, z)Sg{y, z)dx dy 
. Ri R2 

= HI lM^,y))5g(x,yydxdy+^-^ jj Sg(x, z^dx dz + ^-^^^^ JJ 5g(x,zy 



dx dz 



[0,1]2 



ill 



R2 



(30) 



+3A2 [' e [h{zf + f2{zf)] + 2(2e - e)h{z)f2{z) dz 
Jo 
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Note that by Cauchy-Schwarz, 



(31) 



(32) 



1/2 



{5g{x,z)fdx > 2 



1/2 



5g{x,z)dx\ =2h{zf 



1/2 



{5g{x,z)fdx > 2 / 5g{x,z)dx\ =2h{zf. 



1/2 



Since /q (e) and Iq^Ic — e) are positive, 51 is bounded from below by 



/"(e) 



+ 



[0,1]2 

/J(2e - 6) 



1/2 



/i(^)^rf^+ / h{zYdz 

'l/2 



1/2 



/2(^)^rf^+ / h{zYdz 
'1/2 



(33) + 3A2 / rfz[e(/i(z)2 + /2(z)2) + 2(2e-e)/i(z)/2(z)] 
Collecting terms and applying equation ( !26l) . this bound becomes 



1/2 



[0,1]2 



(34) 

where 
(35) 

(36) 

(37) 



I'^{g{x,y))5g{x,yYdxdy + / dz[cih{zY + C2f2{zy + 2c^h{z)h{z)] 

Jo 

+ I dzlcMzY + C2fi{zY + 2cMz)Mz)], 

Jl/2 



Cl 
C2 
C3 



TO 

2 

/^(2e 



e(/^(2e 



2(e-e)2 
_e) ^e(/^(2e-e) 



2 ■ 2(e-e)2 
(2e-6)(/-(2e-6)-/-(6)) 
2fe-e)2 



Note that all coefficients are positive, and that C2 > 1. As e — )■ 0, Ci goes to +oo as 1/e, 
while C3 only diverges as — ln(e). Since C1C2 > C3 for small e, the integrand for each z is 
positive semi-definite, so the integral over z is non-negative, and we obtain 



(3J 



where we used the fact that Iq{u) > 2 for all u. 



□ 



Any global minimizer must be 0(e) close to go, and hence 0(e) close to our specified 
perturbative minimizer. This means that the only way for them to differ is through a 
complicated bifurcation of minimizers at go, despite the uniform bounds on 61 as we approach 
the boundary. The difference between these hypothetical minimizers and go would not be 
pointwise small, but would merely be small in an sense. 
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For instance, consider graphons of the form 

{p X < c < y or y < c < x; 
a x,y<c; 
13 x,y>c, 

where c is a parameter that we will vary and p, a and /3 are constants that depend on c. 
For each c sufficiently close to 1/2, it is possible to find a graphon of this form such that 

g{x,y)g{y, z)g{x, z)dx dy dz = t and / / g{x,y)dx dy = e, and such that the Euler- 




Lagrange equations ( l25l) are satisifed. Call this graphon gc{x,y). A lengthly calculation 
shows that 

(40) "''^'^ 



> 166^ 



c=l/2 

for small t, indicating that (nearly) bipartite graphs with pieces of unequal size have a higher 
rate function than g. This provides strong evidence that our perturbative solution is in fact 
a global minimizer for sufficiently small t. 

Corollary 5.2. Assuming our perturbative solution is the global optimizer, there is a phase 
transition near the boundary point (1/2, 0) between the first and second scallop. 

Proof. Our perturbative solution yields a formula for the entropy: 

(41) 5e, = -^[/o(e) + /o(2e-e)]. 

This formula for the entropy cannot be extended analytically beyond e = (1 + e)/2, as 
d'^s/de^ diverges as e — )■ (1 + e)/2. However, e = (1 + e)/2 corresponds to t = (e^ + 3e)/4, 
or, using the more basic variable e. 



(42) t= [(2e-l)=^ + 3(2e-l)]/4, 

which is in the interior of (e, t) space. (Since the graphon g{x, y) is nowhere zero, it differs in 
form from the graphons describing graphs with minimal t.) Thus Se,t must fail to be analytic 
in some neighborhood of the first scallop. □ 

Of course there must also be a phase transition, presumably from this bipartite phase to 
a homogeneous phase, if one fixes e and raises t, which we see as follows. 

Corollary 5.3. Assuming our perturbative solution is the global optimizer, there is a phase 
transition as one raises t, for any fixed < e < 1/2. 

Proof. Recall from Corollary 15.21 the connection between t and e: 



(43) t = e3 - (e - ef . 

Note that t is an increasing function of e and reaches the value 2e^ when e = 2e. From 
equation fj^Tl) for the entropy we see that it cannot be extended analytically to t > 2e^, yet 
for e < 1/2 we have 2e^ < < e^^"^ so (e, t) is in the interior of the phase space (Figure 
3). □ 
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6. Conclusion 

Our goal is a theoretical understanding of why materials develop a solid phase at high 
enough particle density (pressure) and low enough energy density (temperature). As dis- 
cussed in the introduction, this is presumably 'due to' the repulsive part of the intermolecular 
interaction. Simplifying the Lennard- Jones potential, the interaction V{r) = r~^^ for sep- 
aration r should give the desired effect. Working in the microcanonical ensemble the main 
problem would be to show why, at high particle density, the energy ground state is crys- 
talline and that this survives at low, not necessarily minimum, energy density. This has been 
a famous unsolved problem for many years ^ |Brl lUht 151] lARlj . 

We have taken an unusual path here in order to make some progress on this classic problem, 
working in a model in which there are energy ground states of different type as one varies 
the particle density. Such a model might have a phase diagram as in Figure 4. (A simple 
model using short range forces showing different crystal structures in the ground state in one 
spatial dimension is the 'shift model' [NR] . Of course the crystal structures do not survive 
to higher energy density in one dimensional models such as that.) 



Temperature 




Pressure 

Figure 4. A schematic pressure/temperature phase diagram with different 
crystalline phases 

Indeed, many real materials display a range of crystal structures. Most have a close packed 
crystal structure at high particle density, but many are looser packed at lower density, for 
instance simple cubic for table salt. This feature of varying crystal structure is useful for us 
because there is then a phase transition between the different crystal structures near energy 
minimization, the lower edge of the phase space in Figure 4. This is significant because one 
can use this to analyze the structure of the low energy phases, separately from analyzing 
whether those structured phases melt to a fluid. 

This is the path we have taken in our mean field model. We have shown that our 
(2-parameter, bipartite) graphons g of Theorem 15. II maximize the entropy density at least to 
second order in perturbation theory, among states with given, and low, particle and energy 
densities. Assuming the g are actually global maximizers we then proved the entropy density 
would have to lose analyticity as the density of the state approaches the tripartite regime. 
We also show that the entropy density must suffer a phase transition as the energy density 
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is raised sufficiently high, presumably from the structured bipartite phase to a homogeneous 
high energy phase. 

We expect that a more complicated analysis could produce appropriate graphons g^'^^ k > 
1, near each of the higher density (multipartite) ground states, with a transition near each 
scallop intersection. Intuitively this suggests a mechanism whereby as particle density is 
increased, near the energy ground state, the system progressively transitions through finer 
and finer structure; for high particle density most configurations would consist of many 
interacting 'parts', in a crude approximation to how the unit cells in a crystal behave. 

Our results on phase transitions require that the graphons of Theorem 15. II be in fact global, 
not just local, maximizers of the entropy density. In a future paper we use a symmetry to 
prove that these graphons are indeed the unique global maximizers at least for energy density 
in the range < t < 1/8 and particle density e = 1/2, and we can then see the transition on 
this curve. However we still cannot prove the graphons are the global optimizers of entropy 
density for (e, t) in a two-dimensional region, as is needed to fully justify the notion of a 
structured phase. (See |AR2j for a variant of this approach.) 

In conclusion we emphasize that our key tool was Theorem 13. 11 an optimization formula for 
the asymptotic entropy density, based on the graph limit formalism. The graphon formalism 
is a powerful tool for dealing with the infinite volume limit in mean field models, which we 
have used to make some small progress on the classic problem of the fluid/solid transition. 

Acknowledgements: We gratefully acknowledge useful discussions with Francesco Maggi 
and Peter Winkler. 
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